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Introduction 

Attitude oscillations have been observed during the terminal part 
of automatic landings. A manual landing by the astronaut circumvents 
the attitude oscillations in P64 and P65. But it is desirable to build an 
automatic system from which the crew takes over because of preference 
not necessity. Allan Klumpp and Bill Widnall have concluded that the 
oscillations are due to the lags between the state vector times and the 
realization of the guidance commands by the two control systems, the 
engine throttle servo and the digital autopilot. Allan's engineering 
simulation appears to confirm this. This memo gives a set (which 
ought to turn out to be somehow equivalent to Allan's) of lunar landing 
guidance equations which provide lead compensation for the lags. 

Definition 

Here are some definitions necessary to understand the rest of 

state vector time, the time at which the PIPA's are 
read. 

the position and velocity at t , i. e. r(t Q ), v(t Q ). 
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t = t Q + r, a time r seconds later than state vector time, 

but never greater than T. r is the amount of time the 
guidance equations are projected into the future from 
state vector time. 

T = the time at which the desired position, velocity, and 

acceleration vectors are to be achieved. 
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the position, velocity and thrust acceleration vectors, 
desired at T. 


the n time-to-go". The time between t Q and T. Thus, 



the "diminished" time-to-go, which tries to take into 
account the lags between t Q and the realization of 
guidance commands. 



the gravitation acceleration at state vector time (output 
of AVERAGEG) and at T, respectively. 



the guidance equations commanded thrust acceleration 
r seconds after state vector time. 


three guidance parameters which are functions of T^ q , 
the state at t Q , the desired state at T, and the desired 
thrust acceleration vector at T. These are such that 
the exact solution to the guidance problem is given by 

— TC^ t) = - + (T_t) - + (T ' t>2 - ■ £ (t) 

This equation is valid for the entire future and the 
guidance parameters change only because of navigation 
errors, terrain irregularities, control errors, re- 
targetting. Realistically, the guidance parameters 
change very slowly. 



The Equations 


The equations are given here for those who are interested in the 
computational aspects without wading through the derivation. A second 
memo will give my derivation of the equation. The equations are given 
in steps which the computer program might take. 


1. Determine r^, v q and t Q . That is, read the computer clocks, 
the PIPA's, and perform the AVERAGEG equations. 

2. Compute a, the first guidance parameter 

a = a TD + £(T) 

(Note: This does not have to be computed every two seconds; 
but only when the targetting changes. ) 


3. Compute 


T = T - t 
go o 


4. Compute the following functions of state at T and t^ 

x = v„ - v - T a 
— — D — o go - 

y = r - r - T v - (T 2 / 2)a 

5. Compute the 2nd and 3rd guidance parameters. 

b = (18/T go 2 ) x- (24/T go 3 ) Z 

c - -(24/T go 3 )x+(36/T go 4 ) Z 

6. Compute the command thrust acceleration vector corresponding 


* 

to t . 


T = T - r (if T > 0; otherwise T = 0) 
go go go go 


(t") = a + T"' b + (T* ) 2 c - g 
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(Note: above should really be g(t ); but the g(t) vector changes 

negligibly in 3 seconds, the approximate proposed value 


for r - t - t . ) 

o 



7. 


Compute the desired thrust acceleration magnitude and 
direction. 

Ox 

direction = unit [£ TC (Oj 

magnitude = abval f a TC (t V )J 
and output these commands. 


A few words are in order about the adavntages of these equations. 

When T becomes small the navigation noise on r and v will 
S' - * — o — o 

begin to make b and c go "wild" (because of the rapidly increasing gains 

on x and y in step 5). But a, jb, and c^ should be constants (if navigation 
and control are perfect) and step 6 can be performed with b and c com- 
puted five or ten seconds or even longer in the past. Thus, when T 

go 

becomes less than, say, ten seconds, steps 4 and 5 can be skipped 
and step 6 used to project the desired thrust acceleration vector into 
the future. 

Another consideration which this form of the equations makes 
possible is the relaxation of final position constraint at a certain T 

go 

but the continued maintenance of the velocity constraint at a smaller 

T . 

g° 

Ox 

In steps 6 and 7, I showed the same t for projecting both the 
direction (input to DAP) and magnitude (input to engine throttle servo); 
but for very little additional coding a different T* q and a TC (t*) can 
be generated for the commands to the DAP and throttle. This has 
virtue if the time constants (to which r is related) are very different. 

Note, finally, that we could generate an analytically perfect w^ 
and d (magnitude)/dt. The w D could drive the DAP. 

The above equations make no approximations except for the slight 
one noted in step 6. But I may have blundered algebraically in the 
derivation. The derivation will follow in the next memo and confirm 
or correct the equations. 



